Entanglement Sharing in the Two-Atom Tavis-Cummings Model 
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Individual members of an ensemble of identical systems coupled to a common probe can become 
entangled with one another, even when they do not interact directly. We investigate how this type of 
multipartite entanglement is generated in the context of a system consisting of two two-level atoms 
resonantly coupled to a single mode of the electromagnetic field. The dynamical evolution is studied 
in terms of the entanglements in the different bipartite partitions of the system, as quantified by 
the I-tangle. We also propose a generalization of the so-called residual tangle that quantifies the 
inherent three-body correlations in our tripartite system. This enables us to completely characterize 
the phenomenon of entanglement sharing in the case of the two-atom Tavis-Cummings model, a 
system of both theoretical and experimental interest. 



I. INTRODUCTION 

The control of quantum systems through active mea- 
surement and feedback has been developing at a rapid 
pace. In a typical scenario, a single atom is monitored 
indirectly through its coupling to a traveling probe such 
as a laser beam. The scattered beam and the system be- 
come correlated, and a subsequent measurement of the 
probe leads to backaction on the system. A coherent 
drive applied to the system can then be made conditional 
on the measurement record, leading to a closed-loop con- 
trol model 0] . Such a protocol has been implemented 
to control a single mode electromagnetic field in a cavity 
0], and has been envisioned for controlling a variety of 
systems such as the state of a quantum dot in a solid 0] , 
the state of an atom coupled to a cavity mode 0], and 
the motion of a micro-mechanical resonator coupled to a 
Cooper pair charge box [6j. 

A common theme in the examples given above is that 
measurements are made on single copies of the quantum 
system of interest. However, in many situations one does 
not have access to an individually addressable system. 
In a gas, for example, preparing and/or addressing indi- 
vidual atoms is extremely difficult. In situations such as 
this, it is useful to think of the entire ensemble as a sin- 
gle many-body system. Indeed, recent experiments 00] 
and theoretical proposals 0] have explored the control of 
such ensembles from the point of view of the Dicke model 
(To| . where a collection of N two-level atoms is treated 
as a pseudo-spin with J = N/2. 

Measurement backaction on the pseudo-spin can lead 
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to squeezing of the quantum fluctuations 01 , which 
may be enhanced through active closed-loop control 
0, 0| • This squeezing can reduce the quantum fluctu- 
ations of an observable as in, for example, the reduc- 
tion of "projection noise" leading to enhanced precision 
measurements in an atomic clock [Tlj |. Moreover, spin- 
squeezing is related to quantum ent ang lement between 
the atomic members of the ensemble [l2j . This entangle- 
ment arises not through direct interaction between the 
atoms, but through their coupling to a common "quan- 
tum bus" in the form of an applied probe. 

Measures of entanglement associated with these spin- 
squeezed states have been studied by Stockton, et. al, 
|l0| under the assumption that all of the atoms in the en- 
semble are symmetrically coupled to the bus. However, 
completely quantifying entanglement in the most general 
cases is extremely difficult, and as yet, an unsolved prob- 
lem . In this article we consider the simplest possible 
ensemble consisting of two two-level atoms. Although at 
first sight this might appear trivial, when such a system 
is coupled to a quantum bus a rich structure emerges. 
Again, we consider the simplest realization of the bus 
- a single mode quantized electromagnetic field. The 
resulting physical system then corresponds to the two- 
atom Tavis-Cummings model (TCM) [l4|. A thorough 
understanding of the dynamical evolution of the TCM 
has obvious implications for the performance of quan- 
tum information processing 0, IwL Il7j . as well as for 
our understanding of fundamental quantum mechanics 
[TM Il8[ . Bipartite entanglement has been investigated in 
this system for the one-atom case, known as the Jaynes- 
Cummings model, for initial pure states |l9j | and mixed 
states pa E| of the field. 

Taken as a whole, the two-atom TCM in an overall 
pure state constitutes a tripartite quantum system in a 
Hilbert space with tensor product structure 2 <S> 2 ® oo. 
Entanglement in tripartite systems has been studied by 
Coffman, et. al, [22] for the case of three qubits. They 
found that such quantum correlations cannot be arbitrar- 
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ily distributed amongst the subsystems; the existence of 
three-body correlations constrains the distribution of the 
bipartite entanglement which remains after tracing over 
any one of the qubits. For example, in a GHZ-state, 
| GHZ) = 1/V2(|000) + |111», tracing over one qubit re- 
sults in a maximally mixed state containing no entangle- 
ment between the remaining two qubits. In contrast, for 
aW-state, |W) = 1/V§(|001) + |010) + | 100)), the average 
remaining bipartite entanglement is maximal j^. Coff- 
man, et. al., analyzed this phenomenon of entanglement 
sharing [22^ . using an entanglement monotone known as 
the tangle, a simple generalization of the more familiar 
concurrence 0> HH • They also introduced a new quan- 
tity known as the residual tangle, in order to quantify the 
irreducible tripartite correlations in a three qubit system 

M 

In this article we extend the analysis of entanglement 
sharing to the case of the two-atom TCM. This has im- 
plications for the study of quantum control of ensembles. 
For example, if we imagine that the quantum bus is mea- 
sured, e.g., the field leaks out of the cavity and is then 
detected, then the degree of correlation between the field 
and one of the atoms determines the degree of backaction 
on one atom. We can then quantify the degree to which 
one can perform quantum control on a single member of 
an ensemble even when one couples only collectively to 
the entire ensemble. We will accomplish this by extend- 
ing the residual tangle formalism of Coffman, et. al, to 
our 2 ® 2 ® oo system. 

The remainder of this article is organized as follows. 
The important features of the TCM are reviewed in Sec. 
UTI and the applicable measures of entanglement are in- 
troduced in Sec. II I II With this formalism in hand, we 
calculate the tangle for each of the bipartite partitions 
of this tripartite system in Sec. IIVI We will find an 
approximate analytic expression for the tangle between 
the field and the ensemble in the limit of large average 
photon number and in the Markoff approximation which 
provides further insight into these results. In Sec. we 
study the irreducible tripartite correlations in the system 
using our proposed generalization of the residual tangle. 
Finally, we summarize our results in section FVII 
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where and display a local SU(2) algebra for 
the i-th atom in the two-dimensional subspace spanned 
by the ground and excited states {!<?), |e)}, and a (a') 
are bosonic annihilation (creation) operators for the 
monochromatic field. The Hilbert space TL of the joint 
system is given by the tensor product Ha 1 ® TLa 2 ®7~Lf 
where TLa ± 0~La 2 ) denotes the Hilbert space of atom one 
(two) and TCf is the Hilbert space of the electromagnetic 
field. 

The total number of excitations K — a^a + ^(eri + 

(2) 

o-r+2) is a conserved quantity which allows one to split 
the Hilbert space TL into a direct sum of subspaces, i.e., 
TL = Y^k=q®^k, with each subspace Qk spanned by 
the eigenvectors {|ee, k — 2), \eg, k — 1), \ge, k— 1), \gg, k)} 
of K with eigenvalue k. The analytic form for the time 
evolution operator within a subspace £Ik is given in Zeng, 
et. al., H3- 

It is assumed throughout that the initial state of the 
TCM system is pure. Furthermore, we consider only the 
effects of the unitary evolution generated by Eq. QJ, i.e., 
we do not include the effects of measurement, nor of mix- 
ing due to environment-induced decoherence |28| , so that 
our system remains in an overall pure state at all times. 
Finally, by assuming an identical coupling constant g be- 
tween each of the atoms and the field, the Hamiltonian 
is symmetric under atom-exchange. This invariance un- 
der the permutation group was used by Stockton, et. al., 
|10| to analyze the entanglement properties of very large 
ensembles. We will also make use of this fact in order 
to reduce the number of different partitioning schemes 
that one needs to consider when studying entanglement 
sharing in the two-atom TCM. 



III. THE TANGLE FORMALISM 



II. THE TAVIS-CUMMINGS MODEL 

The Tavis-Cummings model (TCM) [l4| (also known 
as the "Dicke model" [261] ~) describes the simplest funda- 
mental interaction between a single mode of the quan- 
tized electromagnetic field and a collection of N atoms 
under the usual two-level and rotating wave approxima- 
tions. The two-atom (N = 2) TCM is governed by the 
Hamiltonian 



The tangle between two qubits in an arbitrary state 
is defined in terms of the concurrence |24l |25| . For a 
ip) of two qubits, the concurrence is given 



where 



ip ) represents the 'spin-flip' 



pure state 
by C (VO = 

of \ip), i.e., tpj = <7 y <£> a y \ip*), and the star denotes 
complex conjugation in the standard basis. 

The generalization of the concurrence to a mixed state 
p of two qubits is defined as the infimum of the av- 
erage concurrence over all possible pure state ensem- 
ble decompositions of p, defined as convex combina- 
tions of pure states Si — {pi,ipi}, such that p = 
EiP#i)(V>i|- In thi s way, C (p) = in£ Si J2iP* c = 
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infSiJ2iPi yPA^iJ ■ Wootters succeeded in deriving an 
analytic solution to this difficult minimization procedure 
in terms of the eigenvalues of the nonHermitian oper- 
ator pp, where the tilde again denotes the spin-flip of 
the quantum state. The closed form solution for the 
concurrence of a mixed state of two qubits is given by 
C (p) — max{0, Ai — A2 — A3 — A4}, where the Aj's are 
ordered in decreasing order [25|. Rungta, et. ai, ex- 
tended this formalism by introducing an analytic form for 
the concurrence of a bipartite system AB, with arbitrary 
dimensions Da and Db, in an overall pure state |29| . The 
analytic form of this quan tity, dubbed the I-co ncurrence, 
is given by C {iPab) = ^/2v a vb [1 - tr {p 2 A % where p A 
is the marginal density operator obtained by tracing the 
joint pure state over subsystem B, and va and vb are 
arbitrary scale factors. 

Coffman, et. ai, defined the tangle T2 for a system of 
two qubits as the square of the concurrence p^ |. i.e., 



t 2 (p) = max{0, Ai - A 2 - A3 - A 4 }^ 



(2) 



Indeed, we may extend this definition directly to the re- 
sult of Rungta, et. ai, in order to obtain an analytic form 
for the tangle r of a bipartite system in a pure state with 
arbitrary subsystem dimensions, 

r {iPab) = C 2 (tJjab) = 2v A v B [1 - tr (p 2 A )] . (3) 

However, when extending this definition to apply to a bi- 
partite mixed state pab with arbitrary subsystem dimen- 
sions, one must find the infimum of the average squared 
pure state concurrence [3(J 



borne derived an analytic form for r {pab) in the case 
where the rank of pab is no greater than two, 

r {pab) = tr {pabPab) + 2\[i B J [l - tr {p 2 AB )] , (5) 

where pab now represents the universal inversion p9| 
of pab, and A^„ is the smallest eigenvalue of the M 
matrix defined by Osborne [3j . The important point is 
that Eq. JSJ) yields a closed form which, as we will see in 
Sect. IIV CI is directly applicable to a specific bipartite 
partition of the two-atom TCM. 



IV. BIPARTITE TANGLES IN THE TWO-ATOM 

TCM 

Let the two atoms in the ensemble be denoted by A\ 
and A2, respectively, and the field, or quantum bus, by 
F. Because of the assumed exchange symmetry, there 
are four nonequivalent partitions of the two-atom TCM 
into tensor products of bipartite subsystems: (i) the field 
times the two-atom ensemble, F® {A1A2), (ii) one atom 
times the remaining atom and the field, A\ ® {A2F) = 
A2 ® {AiF), (iii) the two atoms taken separately, having 
traced over the field, A\ ® A 2 , and (iv) one of the atoms 
times the field, having traced over the other atom, A\ ® 
F = A2® F '. We calculate how the tangle for each of 
these partitions evolves as a function of time under TCM 
Hamiltonian evolution using the formalism reviewed in 
Sec. |n] We take the initial state to be a pure product 
state of the field with the atoms. We will capture the 
key features of the tangle evolution by considering three 
classes of initial state vectors, 



t{ P ab) = inf ]>>C 2 (V^) 

i 

= 2v A v B v^Y iPi l\-tr (p^Y }, (4) 

where we have used Eq. <(3J for the pure state tangle with 
p^ A as the marginal state for the i th term in the ensemble 
decomposition. 

At this point we note that the scale factors va and vb 
in Eqs. and l@J), which may in general depend on the 
dimensions of the subsystems Da and Db respectively, 
are usually set to one so that agreement with the two 
qubit case is maintained, and so that the addition of ex- 
tra unused Hilbert space dimensions has no effect on the 
value of the concurrence [2^. We will find in Sect. [V] 
when we attempt our own further generalization of the 
tangle formalism, that it is useful to take advantage of 
this scale freedom. For now, however, we adopt the usual 
convention both for the sake of clarity, and to demon- 
strate exactly where in our proposed generalization this 
freedom is required. 

Using the definition given by Eq. (0J), dubbed the I- 
tangle in reference to the work of Rungta, et. al., Os- 



|ee, n) 



|ee, a) or \gg, a) 



(6a) 



(6b) 



and 

(N) + \ge)) ® \a) or -±= {\gg) + \ee)) \a) , (6c) 

where \g{e)) denotes the ground (excited) state of the 
atom, \n) denotes a Fock state field with n photons, and 
I a) denotes a coherent state field with an average number 
of photons given by (n). The alternatives in Eqs. (|6bjl 
and (|6*c|) arise from the fact that, in the limit of large 
(n), the evolution of all of the tangles are found to be 
identical for the two different initial atomic conditions, 
as we will see below. 



A. Field-Ensemble and One Atom- Remainder 
Tangles 

Under the assumption that the system is in an over- 
all pure state, we may easily calculate the tangles in 
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FIG. 1: (Color online) TCM evolution for both atoms ini- 
tially in the excited state and the field in an initial Fock state 
with n = 10. (a) Solid curve (red): Field-ensemble tangle 
t f(a 1 a 2 )\ Large-dotted curve (green): One atom-remainder 
tangle ta 1 (a 2 f)', Dashed curve (black): Atom-atom tangle 
t Ai a 2 (Note that the atom-atom tangle is always zero for 
the given initial condition); Dot-dashed curve (pink): Single 
atom-field tangle ta x f', Dotted curve (blue): Residual tangle 
ta 1 a 2 f- (b) Atomic inversion of the ensemble. 
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FIG. 2: (Color online) TCM evolution for both atoms initially 
in a stretched state and the field in an initial coherent state 
with (n) — 100. (a) Solid curve (red): Field-ensemble tangle 
t f(a 1 a 2 )\ Large-dotted curve (green): One atom-remainder 
tangle ta 1 (a 2 f)\ Dashed curve (black): Atom-atom tangle 
ta 1 a 2 \ Dot-dashed curve (pink): Single atom- field tangle 
ta x f', Dotted curve (blue): Residual tangle ta 1 a 2 f- (b) 
Atomic inversion of the ensemble. 



partitions (i) and (ii) above by applying Eq. J3J), with 
va = vb = 1' Specifically, 

r F(AlA2) = 2 [1 - tr ( P 2 F )] = 2 [1 - tr (p 2 AlA2 )] , (7) 
and 

t Ai (a 2 f) = 2 [1 - tr (p 2 Ai )] = 2 [1 - tr (p^)] , (8) 

where we have used the fact that the (nonzero) eigen- 
value spectra of the two marginal density operators for 



a bipartite division of a pure state are identical | l5l l31| 
in obtaining the rightmost equalities. These tangles have 
implications for the quantum control of atomic ensem- 
bles. Because the overall system is pure, any correlation 
between the field and the ensemble is necessarily in one- 
to-one correspondence with the amount of entanglement 
between these two subsystems. The quantum backaction 
on the ensemble due to measurement of the field is thus 
quantified by Eq. Q. Alternatively, a measurement of 
one atom leads to backaction on the remaining subsystem 
as described by Eq. JSJ). 

The time evolutions for each of the different tangles, 
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FIG. 3: (Color online) TCM evolution for the atoms ini- 
tially in one of the symmetric states and the field in an ini- 
tial coherent state with (n) = fOO. (a) Solid curve (red): 
Field-ensemble tangle t F (a 1 a 2 )', Large-dotted curve (green): 
One atom-remainder tangle ta x (a-,f)', Dashed curve (black): 
Atom-atom tangle ta 1 a 2 ', Dot-dashed curve (pink): Single 
atom-field tangle ta x f', Dotted curve (blue): Residual tangle 
ta 1 a 2 f- (b) Atomic inversion of the ensemble. 

corresponding to the initial conditions given by Eqs. I|tja|) 
- (IrJcl) . are shown in Figs. Qfa) - EI a) respectively. Figs. 
mt>) - EHb) show the time evolution of the atomic inver- 
sion, defined as the probability of finding both atoms in 
the excited state minus the probability of finding both 
atoms in the ground state, for reference purposes. We 
find, under certain conditions, that the two stretched 
states in Eq. (|6bl) lead to identical evolution for the tan- 
gles in all of the bipartite partitions of the system, cor- 
responding to the curves shown in Fig. HJa). Similarly, 
the two symmetric states given in Eq. (jfjc|) both yield the 
curves shown in Fig. Ofla). This behavior can be derived 
under a set of highly accurate approximations. In the 



limit of large average photon number, an initial coher- 
ent state field with zero phase will remain approximately 
separable from the atomic ensemble in an eigenstate of 
J x = J + + J_ up to times on the order of (n) / g where, 
in the pseudospin picture, J± = +a± . This follows 
immediately from the time evolution operator generated 
by Hi n t in Eq. in the interaction picture. The key 
observation is that, for a macroscopic field, the removal 
or addition of a single photon has a negligible effect. This 
allows one to approximate the time evolution operator by 
CX p- lH intt/h exp-^v 7 ^")- 7 *'. Thus, the eigenstates of 
J x form a convenient basis to use in describing the state 
of the atomic ensemble. This approach was taken by 
Gca-Banacloche in analyzing the behavior of the single 
atom Jaynes-Cummings model |32j and the generation of 
macroscopic superposition states |33| . and extended to 
the multi-atom TCM by Chumakov, el. al., [31111 HI! • 
We take as the appropriate basis the three symmetric 
eigenstates of J x , which we label by m — -1, 0, and 1; 
the singlet state, J = 0, is a dark state and thus does 
not couple to the field. Writing the initial state of the 
system as 

l 

iv(o)) = Yl d ™ i m > ® i«) » ( 9 ) 

m— — 1 

and using the factorization approximation |34j . we find 
that the state of the system up to times on the order of 
(n) I g is given by 

l 

\i/j(t))* d m \A m {t))®\4> m {t)), (10) 

m=—X 

where |A m (t)) and \4> m (t)) are the time-evolved atomic 
and field states, respectively. The marginal density op- 
erator for the two atoms is then 

PMA 2 (t) « (Mt)\ fmi (gt, (n)) , (11) 

I, in 

where f m i (gt, (n» = (n | <j) m (t)) (cj>i(t) \n). We find 
that this function has "memory" only for t <C yj (n) /g, 
and behaves very much like a delta function for longer 
time scales. Effectively, the large dimensional Hilbert 
space of the field acts as a broadband reservoir for the 
atoms - the generalization of the familiar "collapse" 
phenomenon in the Jaynes-Cummings model. This 
"Markoff" approximation is valid up to times on the or- 
der of 2nyJ (n) I g, corresponding to the well-known re- 
vival time in the Jaynes-Cummings Model [22|. Making 
this approximation in Eq. (|11|) . the states \A m ( t)) act 
effectively as a "pointer basis" for decoherence j2a| of the 
atomic density matrix, i.e., 

PMA 2 (t) w Y l d ™! 2 I^W) ^M*)l • (12) 

m 
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Substituting this formula into Eq. J7J yields 

T F(AlA2) {t)^ 2^1 -\[c-h(t')]Y (13) 

where 

c = 4(|d_ 1 | 4 + Mo| 4 + Mi| 4 ) +2 |d | 2 Mil 2 

+ |d_ 1 | 2 (2|d | 2 + 3|d 1 | 2 ) -4|^ 1 | 2 |d 1 | 2 (14) 

h(t') = 2 |d | 2 (M-i| 2 + Mil 2 ) cos(40 

+ |d_i| 2 |di| 2 cos(80, (15) 

and 

f= , gt (16) 
V<»> - T + I 

Under the factorization and Markoff approximations, 
the field-ensemble tangle is given by a constant term c 
that depends only on the initial probabilities to find the 
atomic ensemble in each of the J x eigenstates, and a time- 
dependent piece h(t'). These probabilities depend solely 
on the absolute squares of the expansion coefficients of 
the initial atomic state given by Eq. ©. It is now clear 
why certain initial atomic conditions result in identical 
evolution for the different tangles. For example, the ini- 
tial atomic states \gg) and |ee) both satisfy 

\d-i\ = |di| = ~ and |d | = (17) 

corresponding to identical evolution for all of the tangles 
shown in Fig. Ufa). Similarly, the initial atomic states 
1/V2{\eg) + \ge)) and 1/V2(|ee) + \gg)) both satisfy 
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FIG. 4: (Color online) Exact field-ensemble tangle: Solid 
(red) curve, and approximate formula: Dashed (black) curve, 
as functions of time. 

is seen to track the exact evolution extremely well over 
the range of its validity. The discrepancy at very small 
times is explained by the fact that at these times the 
Markoff approximation breaks down. It is also seen that 
our approximate solution does not capture the small dip 
in the field-ensemble tangle occurring at t = H\J (n) / g. 
The absence of this feature can be explained by noting 
that in making the Markoff approximation we have ef- 
fectively wiped out any information regarding the initial 
coherence between the \m = —1) and \m = 1) states. The 
presence of this dip is then seen to be dependent upon 
the existence of this coherence. This is borne out by the 
fact that the dip in the field-ensemble tangle in Fig. |2Ja) 
is much shallower than that in Fig. Ofa), where the ini- 
tial atomic expansion coefficients are given by Eqs. (|17|) 
and ljl%|) . respectively. 

B. Atom- Atom Tangle 



|d_i| = Mi| = -j^and |do| = 0, (18) 



corresponding to the curves shown in Fig. EJa). More 
generally, this property holds for any class of initial states 

(0)) having the form of Eq. © such that 



m G {—1,0,1}. One immediate consequence of 

this result is that the relative phase information encoded 
in the initial state of the atomic system is irrelevant to 
the evolution of the field-ensemble tangle. 

The field-ensemble tangle calculated according to Eq. 
@ and the approximation given by Eq. (|13f) for an initial 
stretched atomic state and an initial coherent state field 
with (n) = 500 are shown by the solid (red) and dashed 
(black) curves in Fig. respectively. The approximation 



Given an initial state, we time-evolve the system ac- 
cording to the dynamics governed by Eq. l[TJl. and then 
trace over the field subsystem. The tangle of the two- 
atom mixed state p Al A 2 {t) may then be calculated ac- 
cording to Eq. J5J). The resulting atom- atom tangles 
corresponding to the initial conditions in Eqs. (|6a|l - 
(|dc|) are depicted by the dashed (black) curves in Figs. 
M a ) - EI a ) > respectively. These curves yield direct in- 
sight into the state of the atomic ensemble as a function 
of time. Specifically, the atom-atom tangle quantifies the 
degree to which the ensemble behaves as a collective en- 
tity, rather than as two individual particles. 

It is somewhat surprising that for the initial condition 
given by Eq. i|tja|) . *- e -; when the field is initially in a 
Fock state with any value for n, the atom-atom tangle 
remains zero at all times, whereas the evolution of the 
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atom-atom tangle resulting from an initial coherent state 
field is nontrivial and, in general, nonzero. As a first 
step towards understanding these observations, we have 
studied the evolution of the atom-atom tangle for other 
initial conditions. When the field is initially in a Fock 
state and both atoms start in the ground state, the loss 
of an excitation in the field can result in the creation 
of an excitation in the atomic ensemble. This produces 
entanglement between the field and the ensemble and in 
the single atom- field and one atom-remainder partitions. 
Since it is not possible to distinguish in which atom the 
excitation is created, the two atoms become entangled 
with each other as well. It is found that the atom-atom 
entanglement falls off as 1/n 2 so that, in the limit of a 
highly excited Fock state, these initial conditions yield 
results reminiscent of those found in Fig. ^a). Specifi- 
cally, we find that the entanglement in all of the different 
subsystem partitions always oscillate in phase at twice 
the Rabi frequency, and that the atom-atom tangle ap- 
proaches zero as n becomes large. 

Next, we considered the case when both atoms initially 
reside in a stretched state, and the initial field state con- 
sists of a coherent superposition of two neighboring Fock 
states. We find, on a time scale much longer than that 
given by the inverse of the associated Rabi frequencies, 
that the overall behavior again closely resembles the evo- 
lution seen for an initial field consisting of a single Fock 
state. Specifically, we find that the general features of all 
of the different bipartite tangles oscillate in phase with 
one another. However, on much shorter time-scales, the 
effects of dephasing between the two Rabi frequencies be- 
come apparent, yielding the first clues regarding how the 
observed coherent state behavior arises in terms of initial 
Fock state superpositions. For example, it seems likely 
that these observations provide insight into the Fock-like 
behavior seen in Fig. |2Ia) at the revival time, when there 
is a partial rephasing of the Rabi frequencies. At this 
time all of the bipartite tangles decrease simultaneously, 
while at other times the tangles in certain bipartite par- 
titions may be completely out of phase with one another. 
We are currently performing a detailed study of the en- 
tanglement that can be dynamically generated between 
the two atoms under TCM evolution for the most gen- 
eral initial conditions in order to better understand this 
behavior. 



C. Single Atom-Field Tangle 

The final bipartite partition of the two-atom TCM is 
that consisting of a single atom, say A\ , as one subsystem 
and the field F as the second subsystem. Again by ex- 
change symmetry ta 1 f = ta 2 f, so we need calculate only 
one of these quantities. Because the tripartite system is 
in an overall pure state, the Schmidt decomposition the- 
orem |l5l l3l| implies that the marginal density operator 
Pa x f has at most rank two. The rank of the reduced 
density matrix is set by the dimension of the smallest 



subsystem, which in this case is a two-level atom. This 
is exactly the scenario envisioned by Osborne [3(| , as de- 
scribed in Sec. IIIII The tangle corresponding to this 
partition, A\ ® F, is computed by first tracing over the 
state of the remaining atom, A2, and then applying Eq. 

Employing this procedure, 

r Al F = tr (p Al FPA lF ) + 2A^„ F) [l - tr (p 2 AiF )] , (19) 

where A„^.„ represents the minimum eigenvalue of the 
Osborne M matrix 30] generated from the marginal den- 
sity operator pa x f- The dot-dashed (pink) curves in Figs. 
Wi a ) -EUa) give the time evolution of the single atom-field 
tangle for the different initial conditions considered. 

We are now in possession of closed forms for the tan- 
gles of all bipartite partitions of the two-atom TCM. Any 
other entanglement that the system may possess must 
necessarily be in the form of irreducible three-body quan- 
tum correlations. In section^we review the residual tan- 
gle formalism introduced by Coffman, et. al, in order to 
quantify this type of tripartite entanglement in a system 
of three qubits. We then propose a generalization of this 
quantity that is applicable to a 2 ® 2 ® D system in an 
overall pure state. This extension of the tangle formal- 
ism allows us to study the phenomenon of entanglement 
sharing in the two-atom TCM. 

V. ENTANGLEMENT SHARING AND THE 
RESIDUAL TANGLE 

The concept of entanglement sharing studied in [22ll37| 
refers to the fact that entanglement cannot be freely dis- 
tributed among subsystems in a multipartite, i.e., tripar- 
tite or higher, system. Rather, the distribution of entan- 
glement in these systems is subject to certain constraints. 
As a simple example, consider a tripartite system of three 
qubits labeled A, B, and C. Suppose that qubits A and B 
are known to be in a maximally entangled pure state, e.g., 
the singlet state, given by \^ab) = l/y/2 (|01) - 1 10) ) 
when written in the logical basis. In this case, it is ob- 
vious that the overall system ABC is constrained such 
that no entanglement may exist either between A and C 
or between B and C. Otherwise, tracing over subsystem 
C would necessarily result in a mixed marginal density 
operator for AB in contradiction to the known purity of 
the singlet state. 

Coffman, et. al., analyze the phenomenon of entangle- 
ment sharing for a system of three qubits in an overall 
pure state in full generality by introducing a quantity 
known as the residual tangle [22|. This definition is mo- 
tivated by the observation that the tangle of A with B 
plus the tangle of A with C cannot exceed the tangle of 
A with the joint subsystem BC, i.e., 

tab + t~ac < T A (BC)- (20) 
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Here, tab and tac are calculated according to Eq. J3J|, 
and ta(bc) ma y be obtained from Eq. 

The original proof [2^ of the inequality in Eq. I)2()[l. 
which forms the heart of the phenomenon of entangle- 
ment sharing for the case of three qubits, may be sub- 
stantially simplified by making use of certain results due 
to Rungta, et. al. Specifically, we note that [2IJ 

tr (p xy p xy ) = l-tr (pi) - tr (p 2 y ) + tr (p 2 xy ) > (21) 

for subsystems x and y having arbitrary Hilbert space 
dimensions. Under the assumption that x and y are in 
an overall pure state with a third subsystem z, Eq. 121|) 
may be rewritten 

tr (p xy p xy ) = l-tr (pi) - tr (p 2 y ) + tr (p 2 z ) > 0, (22) 

where we have used the equality of the nonzero eigenvalue 
spectra of p xy and p z . Then, by the observation [22] that 
for an arbitrary state of two qubits A and B, the following 
upper bound on the tangle defined by Eq. J2J) holds 

t 2 (pab) < tr (pabPab) , (23) 

and by Eq. 10 with vj± = vb = 1, the inequality in Eq. 
(|2U|I follows immediately. 

Subtracting the terms on the left hand side of Eq. l(2T>|) 
from that on the right hand side yields a positive quantity 
referred to as the residual tangle tabc, i- e -, 

TABC = Ta(BC) - T AB - TAC- (24) 

The residual tangle is interpreted as quantifying the in- 
herent tripartite entanglement present in a system of 
three qubits, i.e., the entanglement that cannot be ac- 
counted for in terms of the various bipartite tangles. This 
interpretation is given further support by the observation 
that the residual tangle is invariant under all possible 
permutations of the subsystem labels • 

We wish to generalize the residual tangle, defined for a 
system of three qubits, to apply to a 2 (g> 2 (g> D quantum 
system in an overall pure state so that we may study 
entanglement sharing in the two-atom TCM. Note that 
we already have all of the other tools needed for such 
an analysis. Specifically, from section IIVI we know the 
analytic forms for all of the different possible bipartite 
tangles in such a system. 

Any proper generalization of the residual tangle must, 
at a minimum, be a positive quantity, and be equal to 
zero if and only if there is no tripartite entanglement in 
the system, i.e., if and only if all of the quantum corre- 
lations can be accounted for using only bipartite tangles. 
It should also reduce to the definition of the residual tan- 
gle in the case of three qubits. Further it is reasonable 
to require, if this is to be a true measure of irreducible 
three-body correlations, that symmetry under permuta- 
tion of the subsystems be preserved, and that it remain 



invariant under local unitary operations. Finally, we con- 
jecture that this quantity satisfies the requirements for 
being an entanglement monotone [33,1111 under the set of 
stochastic local operations and classical communication 
(SLOCCs), or equivalently, under the set of invertible 
local operations [23j . We limit the monotonicity require- 
ment to this restricted set of operations since, in the con- 
text of entanglement sharing, we are only concerned with 
LOCCs that preserve the local ranks of the marginal den- 
sity operators such that all subsystem dimensions remain 
constant. 

Let A and B again be qubits, and let C now be a D- 
dimensional system with the composite system ABC in 
an overall pure state. We note that, under these assump- 
tions, we are still capable of evaluating each of the terms 
on the right hand side of Eq. H24|) analytically using the 
results of section ITVl However, we cannot simply use the 
definition of the residual tangle (with C now understood 
to represent a D-dimensional system) as the proper gen- 
eralization for two reasons. First, since the three subsys- 
tems are no longer of equal dimension, symmetry under 
permutations of the subsystems is lost. However, as we 
will see, this problem is easily fixed by explicitly enforc- 
ing the desired symmetry. The second, and more difficult 
problem to overcome is the fact that the inequality given 
by Eq. (|20|) no longer holds for our generalized system 
because \ m in in Eq. JSJ can be negative, implying that 
Eq. H24f) can also be negative. 

The required permutation symmetry may be restored 
by taking our generalization of the residual tangle, which 
we dub the I-residual tangle in reference to previous work 
and continue to denote by tabc, to be 

TABC = g (TA(BC) + T B(AC) + T C{AB)) 
2 

-g {tab + tac + tbc) (25) 

The definition in Eq. I|25|l is obtained by averaging over 
all possible relabelings of the subsystems in Eq. i|24[l . 
By inspection, it is obvious that Eq. I|25[l preserves per- 
mutation symmetry. However, it still suffers from the 
problem that its value can be negative. In order to deal 
with this difficulty, we make use of the arbitrary scale 
factors appearing in Eqs. © and (0J, and discussed in 
section IIIII 

Let d be the smaller of the two 'dimensions' of two 
arbitrary dimensional subsystems x and y, i.e., d = 
min {D x , D y }. Note that by dimension we do not nec- 
essarily mean the total Hilbert space dimension of the 
physical system under consideration, but only the num- 
ber of different Hilbert space dimensions that contribute 
to the formation of the overall pure state of the system. 
This is a subtle but important point which automatic ally 
enforces insights like those due to Rungta, et. al, |29| 
and Verstraete, et. at, 40] which state that the scale 
chosen for a measure of entanglement must be invariant 
under the addition of extra, but unused, Hilbert space 
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dimensions. The two-atom TCM provides one example 
of the relevant physics underlying these ideas. 

Consider, for example, the bipartite partitioning of 
the TCM into a field subsystem with Dp = oo, and 
an ensemble subsystem consisting of the two qubits with 
Da x a 2 = 4. Any entangled state of the overall system 
has a Schmidt decomposition with at most four terms, 
implying that the field effectively behaves like a four- 
dimensional system. Further, since the Tavis-Cummings 
Hamiltonian given by Eq. Q does not induce couplings 
between the field and the singlet state of the atomic en- 
semble, i.e., the singlet state is a dark state, the field be- 
haves effectively as a three-level system, or qutrit, in the 
context of the TCM. Accordingly, in any entangled state 
of the field with the ensemble, the field is considered to 
have a dimension no greater than three. We employ this 
revised definition of dimension throughout the remainder 
of the paper. 

We now make the choice 

v A v B = d/2, (26) 

when calculating each of the bipartite tangles appear- 
ing on the right hand side of Eq. (|25[l . This choice is 
made for several reasons. First of all, it is in complete 
agreement with the two qubit case, yielding vavb = 1 as 
required. Indeed, when A, B, and C are all qubits, the 
residual tangle given by Eq. I|24[) is recovered. Secondly, 
it takes differences in the Hilbert space dimensions of the 
subsystems into account when setting the relevant scale 
for each tangle. This is important since, in order to study 
the phenomenon of entanglement sharing, the tangles for 
each of the different bipartite partitions must be com- 
pared on a common scale. It is reasonable that this scale 
be a function of the smaller of the two subsystem dimen- 
sions since, for an overall pure state, it is this quantity 
that limits the number of terms in the Schmidt decompo- 
sition. Finally, it is conjectured that a proper rescaling 
of the various tangles will result in the positivity of Eq. 

Note that when applying the proposed rescaling to the 
terms on the right hand side of Eq. 1|25[1 . the only term 
affected is tq(ab), which is rescaled by one-half of the 
smaller of the two subsystem dimensions Dc and Dab- 
Each of the other terms remains unaltered since, in each 
case, at least one of the two subsystems involved is a 
qubit. The net effect of this rescaling is to increase the 
'weight' of the tangle between C and AB relative to that 
of the rest of the tangles. This is reasonable when one 
recognizes that both AB, a system of two qubits, and 
C, a D-dimensional system (in the case D > 2), have 
entanglement capacities |37j | exceeding that of a single 
qubit. 

The requirement that the I-residual tangle be invariant 
under local unitary operations follows trivially, since each 
term on the right hand side of Eq. (|25|l is known to sat- 
isfy this property individually. It is still an open question 
as to whether or not the proposed rescaling is sufficient to 
preserve positivity when generalizing the residual tangle, 



Eq. (|24|l . to the I-residual tangle. Eq. I|25|) . However, nu- 
merical calculations give strong evidence that this is the 
case. The I-residual tangle has been calculated for over 
two-hundred million randomly generated pure states of a 
2 ® 2 (g) 3 system and of a 2 £g> 2 ® 4 system, the only non- 
trivial possibilities. In each instance the resulting quan- 
tity has been positive. We conjecture that the I-residual 
tangle satisfies the requirements of positivity and mono- 
tonicity under SLOCC not only for a 2 ® 2 ® D system, 
where closed forms currently exist for all of the terms on 
the right hand side of Eq. I|25|). but for the most gen- 
eral Da <8> Db <8> Dc dimensional tripartite system in an 
overall pure state (with the proper scaling of each term 
again given by Eq. (|26fl ). The I-residual tangle arising in 
the context of the two-atom TCM is shown by the blue 
curves in Figs, ^a) -Eta). 

The residual tangle, as well as our proposed general- 
ization of this quantity, may be interpreted as the irre- 
ducible tripartite entanglement in a system since it can- 
not be accounted for in terms of any combination of bi- 
partite entanglement measures [2^. A slightly different, 
and possibly more enlightening interpretation is that the 
I-residual tangle quantifies the amount of 'freedom' that 
a system has in satisfying the constraints imposed by the 
phenomenon of entanglement sharing. If the I-residual 
tangle of a tripartite system is zero, then each bipar- 
tite tangle is uniquely determined by the values of all 
of the other bipartite tangles. Alternatively, if tabc is 
strictly greater than zero, then the bipartite tangles en- 
joy a certain latitude in the values that each may assume 
while still satisfying the positivity criterion. The larger 
the value of the I-residual tangle, the more freedom the 
system has in satisfying the entanglement sharing con- 
straints. This reasoning highlights the relationship be- 
tween entanglement sharing and the I-residual tangle. 

Finally, we may interpret the I-residual tangle as the 
average fragility of a tripartite state under the loss of 
a single subsystem. That is, if one of the three sub- 
systems is selected at random and discarded (or traced 
over) , then the I-residual tangle quantifies the amount of 
bipartite entanglement that is lost, on average. It is the 
existence of physically meaningful interpretations such as 
these which prompt us to postulate this new measure of 
tripartite entanglement for a 2 ® 2 ® D system in an over- 
all pure state, rather than to rely on previously defined 
measures based on normal forms |40j or on the method 
of hyperdeterminants |4l|. for example. At this point 
it is unclear what, if any connection these entanglement 
monotones have to the entanglement that exists in dif- 
ferent bipartite partitions of the system, a key ingredient 
in any discussion of entanglement sharing. 

The constraint imposed by entanglement sharing on 
the values of the various bipartite tangles, each of which 
is known to be a positive function, is simply that Eq. (|25|) 
cannot be negative. It then follows that the strongest 
constraint of this form is placed on the two-atom TCM 
when the I-residual tangle is equal to zero. This occurs 
(to a good approximation) periodically in Fig. ^a) for 
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the initial condition given by Eq. I|6a[) . It is at these 
points that each bipartite tangle is uniquely determined 
in terms of the values of all of the other bipartite tangles. 
Conversely, at one half of this period when the I-residual 
tangle achieves its maximum value, the various bipar- 
tite partitions enjoy their greatest freedom with respect 
to how entanglement may be distributed throughout the 
system while still satisfying the entanglement sharing 
constraints. The distribution of correlations is, of course, 
still determined by the initial state of the system and by 
the TCM time evolution, both of which we consider to 
be separate constraints. 

Similarly, the dotted (blue) curves in Figs. Eta) and 
Eta) show the evolution of the residual tangle for the 
initial states given by Eqs. I|6bj) and l|Hc| . respectively. 
Note how the more complicated behavior resulting from 
an initial coherent state field arises from a specific su- 
perposition of Fock states, the tangles of which all have 
a simple oscillatory evolution. This suggests that the 
phenomenon of entanglement sharing may offer a useful 
perspective from which to investigate the way in which 
the coherent state evolution results from a superposition 
of Fock state evolutions. 

The fact that the TCM Hamiltonian leads to a nonzero 
I-residual tangle is interesting in its own right. Inspec- 
tion of Eq. shows that this model does not include a 
physical mechanism, e.g., a dipole-dipole coupling term 
enabling direct interaction between the two atoms in the 
ensemble, but only for coupling between the field and the 
atoms. Consequently, all interactions between the atoms 
are mediated by the electromagnetic field via the ex- 
change of photons, and are in some sense indirect. This, 
however, turns out to be sufficient to allow genuine tri- 
partite correlations to develop in the system as evidenced 
by values of the I-residual tangle that are strictly greater 
than zero. 

Finally, we note that we have considered an alternative 
approach to understanding the constraints on the distri- 
bution of entanglement among the different subsystem 
partitions of the two-atom TCM by using the relative 
entropy of entanglement |38j as our entanglement mea- 
sure. This quantity generalizes in a straightforward man- 
ner to the multipartite case ^3 > an d has a clear physical 
interpretation relating the amount of entanglement in a 
state to its distance from the set of separable states [3^ . 
The existence of upper and lower bounds in the tripar- 
tite case yields another method by which to investi- 
gate the genuine three body entanglement arising in the 
two-atom TCM. The results of the relative entropy of en- 
tanglement approach will be presented in a subsequent 
article. 



VI. SUMMARY AND FUTURE DIRECTIONS 

The two-atom Tavis-Cummings model provides the 
simplest example of a collection of two-level atoms, or 



qubits, sharing a common coupling to the electromag- 
netic field. A detailed understanding of the evolution 
of entanglement in different bipartite partitions of this 
model is valuable for both fundamental theoretical in- 
vestigations, and for accurately describing the behavior 
of certain nontrivial, yet experimentally realizable sys- 
tems. Our proposed generalization of the residual tangle 
augments the current formalism, and allows one to ana- 
lyze the irreducible three-body correlations that arise in 
a broader class of tripartite systems, providing a tool use- 
ful for studying the phenomenon of entanglement sharing 
in the context of a physically relevant and accessible sys- 
tem. 

In future work we hope to generalize this analysis to in- 
clude ensembles with an arbitrary number of atoms. This 
will require further extensions of the tangle formalism in 
order to quantify both the entanglement in a mixed state 
of a bipartite system of arbitrary dimensions having a lo- 
cal rank greater than two, and the multipartite entangle- 
ment in a system composed of more than three subsys- 
tems. Ultimately, we hope to connect this analysis to the 
phenomenon of quantum backaction on individual parti- 
cles when the whole ensemble is measured. The tradeoff 
between the information gained about a system and the 
disturbance caused to that system is certainly fundamen- 
tal to quantum mechanics [ij, l45l |46| . However, the re- 
lationship of this tradeoff to multiparticle entanglement 
is far from clear. Such an understanding would not only 
be a crucial step in designing protocols for the quantum 
control of ensembles, but would also provide deeper in- 
sight into the nature of the correlations at the heart of 
quantum mechanics |47| . 
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